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PREFACE 


Space systems planners — space engineers as well as physicists — require for 
their work accurate values of the fluxes and energy spectra of protons and electrons 
trapped in the geomagnetic field. Numerous measurements have been made of the 
Van Allen belt radiations on many satellites, with different instruments and at dif- 
ferent times. Most of the experimenters organize their data in terms of McHwain’s 
B , L coordinates, but intercomparison of the results from different experiments is 
still complicated by temporal variations of the fluxes and differences in energy re- 
sponse of the various instruments used to obtain the data. A need existed, there- 
fore, for collating the individual measurements and for constructing a model 
environment. 

A program, sponsored jointly by the National Aeronautics and Space Administra- 
tion (Defense Purchase Request W -11, 683) and the United States Air Force (Contract 
AF 04(695)-469) is in progress under the direction of Dr. Vette for the express 
purpose of constructing such a model environment. This report contains a descrip- 
tion of the model which has been developed for the inner radiation belt. The graphs 
in this report give a picture of the radiations; card decks containing the actual flux 
values have also been prepared and can be used in subroutines of larger computer 
programs. In addition, codes have been developed for integrating the flux along 
satellite orbits. Persons on NASA sponsored projects may request such integra- 
tionsfromWilmont N. Hess, Goddard Space Flight Center, and those on Air Force 
projects should request it from James I. Vette, Aerospace Corporation. 

As this effort continues, major emphasis will be placed on describing the outer 
radiation belt, updating the model of the inner belt, and developing the capability to 
predict changes that are likely to occur in the particle populations. Hopefully, such 
changes can be correlated with geophysical indicators or other events. 

This effort would not be possible without the assistance given to this project by 
the experimenters who made the measurements originally. All users of this en- 
vironment will greatly appreciate these efforts. It is to be hoped that this original 
enthusiasm will continue and permit us to maintain an up-to-date model environment. 


A. W. Schardt 

Physics & Astronomy Programs 
NASA Headquarters 
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INTRODUCTION 


For 16 months we have been engaged in a NASA/ Air Force supported program to produce 
trapped radiation model environments that will be of general use to space engineers, systems 
planners, and space physicists. An outline of this effort has been described previously (Reference 
1) and several model environments for the inner zone have been distributed. The present report 
will provide (1) more details concerning the production of these environments and (2) a detailed 
comparison with the existing satellite data. 

It is necessary to combine data from many satellites in order to obtain a good over-all pic- 
ture of the trapped radiation. Furthermore, the inclusion of newer data at regular intervals is 
mandatory because our knowledge of the radiation belts increases as improved instruments are 
flown and new ideas are advanced. The appropriate time scale for updating still seems to us to be 
about six months. However, during the initial phase of the present program a longer period of 
time has been required to collect the data, assess them, and produce an output. The inner zone elec- 
tron map AE1 was completed after eight months and the proton maps API through AP4 required an 
additional eight months. However, various procedures had to be learned, computer programs had 
to be developed, and personnel had to be acquired during this period. It is felt that future environ- 
ments can be produced more rapidly and detailed descriptive reports such as this can follow 
closely the initial distribution of the flux maps. 


GENERAL PROCEDURE 

The three quantities that need to be determined by comparing data from various sources are 
the flux distribution, the energy spectrum, and the absolute value. In general, the flux is a function 
of particle energy, time, position, and magnetic or solar activity. A detailed representation in all 
of these variables would not be very practical because it would be too cumbersome to use. Emphasis 
has been placed on dividing the environment into various segments of space, energy, and time so 
that simplifications of representation and use can be achieved. 
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We have started by analyzing the inner zone radiations first for the following reasons. The 
most intense regions of trapped radiation exist there, many satellite programs are concerned with 
the proper design to function reliably in these radiations, and our knowledge of this region is more 
complete. Over much of the inner zone theB, L coordinate system offers a good way to present the 
data. The general variations are long term decay or growth with characteristic times of months 
and years out to L % 2.0; at larger values definite responses to magnetic activities are observed, 
and changes are produced in much shorter times. 

At present, data are received in reports or private communications in a wide variety of for- 
mats. Omnidirectional measurements above a certain energy or within a given energy band are 
common. More recently, unidirectional fluxes perpendicular to the magnetic field line, j i , have 
been measured by instruments on Injun III, Relay I, and Relay II. Other directional data are given 
in terms of equatorial or local pitch angle distributions. Within this framework of measurements 
the data are presented in many ways and coordinate systems such as: 

1. b-L flux maps 

2. r-\ flux maps 

3. FIux-b maps 

4. Flux-x maps 

5. Flux-B/B 0 maps 

6. Local pitch angle distributions 

7. Equatorial pitch angle distributions 

We have chosen to convert all of these into one common format which can be stored on punched 
cards for further processing and for comparison. For the inner zone data, which we consider 
below L = 3.0, we have found the most useful format to be omnidirectional flux plotted versus B. 
Thus all data received is converted to omnidirectional flux. A series of values of flux, f , and B 
for a given l value are punched on cards; enough points are used so that values at other b points 
can be obtained by interpolation. The L values 

L d = 1.1 + O.li ; (i = 1, 2, ••• 19) 


have been chosen as a compromise between processing too much data and obtaining poor inter- 
polations at intermediate L values. As the techniques become more automated, a finer L grid will 
become feasible. 

Hopefully in the future, data can be provided to us in formats which are easy to handle or will 
be provided on punched cards in our standard format. Such a procedure would result in a consider- 
able saving of time and the lag between receipt of data and production of the environment would be 
minimized. 
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, If spectral measurements are available, then a model spectrum can be devised. On the basis 
of this spectrum, all the measurements are converted to a common energy and a detailed com- 
parison results in the selection of the B dependence of the flux. However, the spectrum may have 
to be obtained by an iteration with selected measurements. Comparison with other data can be 
done to determine the feasibility of the model. After the flux-B curves are drawn, a b-l flux map 
is constructed by cross plotting. This presentation is very useful in showing some of the incon- 
sistencies that developed in selecting the B dependence. An iteration between the B-L map and the 
F-B plot generally produces a smooth representation in both systems that is consistent with the 
data. R-\ maps and equatorial pitch angle distributions are constructed from the resulting model 
by further calculation. Finally the model is read into the computer so that integrations over 
selected circular orbits can be performed in a manner in which the spectral information is preserved. 

In the following sections these general procedures will be described in more detail and the 
model environments will be compared with the data used in making them. 


OMNIDIRECTIONAL - UNIDIRECTIONAL FLUX CONVERSIONS 

Since a large amount of input data was received as perpendicular flux, equatorial or local 
pitch angle distributions, it was necessary to develop an efficient computer program to convert 
these data to omnidirectional fluxes. The inverse transformation is also useful in analyzing data 
and in presenting final results. 

The conversion from unidirectional to omnidirectional flux was formulated in a manner very 
similar to that given recently by Roberts (Reference 2). The omnidirectional flux, j , at a given 
position on the L shell given by the field strength, B , in terms of the unidirectional flux, j , at that 
point is 


J(B> L) 


47 T 



j(B, L, d/x , 


(1) 


where 

fj. = cosine of the local pitch angle 

= cutoff value of n which defines the loss cone, i.e., j(B,L,yz) = 0 for . 

In terms of Robert's variable 

x = (1-B 0 /B) ,/2 . 

where B 0 is the field strength at the magnetic equator, the mirror equation is written as 

(l-* 2 )(l-M 2 ) = (l-x»)(l- Ml a ) . (2) 
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If one relates the general values x , ^ to the equator values by means of (2), then 


1 - M 2 


l ~ y 2 

1-X2 


( 3 ) 


where we have written y as the cosine of the equatorial pitch angle. In terms of x and y, (1) becomes 


T , . = 4?r f yc y j(y) dy 

J( X > - (i- x J)l /2 jy =x ( y 2 _ x 2 )l /2 


(4) 


We have suppressed the L variable in (4) and write j(y) for the equatorial pitch angle distribution. 
Of course y c is the cutoff value defining the loss cone. 

The variables x and y cover the range 


0 < x < y c 
0 < y < y c 

which we divide into intervals by means of the points 


y - o, y r y 2 , •••y i • • • y n 

x = 0, x r x ^ t x, ••• x n 


with x. = y. and y n = y c . For a given value of x - x k _ 1 we can write (4) as 


JK-i) 



y j(y) dy 


If the intervals are taken small enough, j will not vary appreciably over the interval and can 
be removed from the integral. Then one obtains 


JK-i) 



(5) 


The quantity y, is the value of y between y._j and y 1 which is most appropriate to use to evaluate 
j and will be discussed later. The integral in (5) is trivial and can be done immediately. If 
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we define 


jK-.) 


J ; 


J(y.) 


(,--^ i) w [(y, 5 '-%;,) vi -(y, )■'*] ■ 


i > k 


i < k 


we can write (5) as the matrix equation 


h 



i i * 


This matrix is similar to that given in Reference 2; however, we have formulated it for unequal 
intervals because the function j(y) changes rapidly as y -y c in many cases and one needs to take 
small intervals in this region. Although this matrix is independent of L, we have found it necessary 
to make the spacing, or selection of y i and x. , depend on L because the cutoff of the fluxes is 
L-dependent. We have made 19 differently spaced 45 by 45 triangular matrices to use at the L values 
1.2, 1.3 ... . 3.0 where the various data have been tabulated for the inner zone. 

We have determined a value for y . in the following manner. The first few terms in (5) are the 
main ones because the integrand is largest for a given interval size and j is generally a mono- 
tonically decreasing function of y . For those regions where j varies slowly with y , setting y equal 
to any point in the interval will give nearly the same answer. However where j changes rapidly, 
the accuracy of this approximate method depends on the choice of y i . Consider the integral for the 
first interval 


I 



j(y) y dy 

( y 2 _ x 2 ) l /2 


and define y t by the equation 

j(y,) [o 1 ") 2 - * 2 ]' /2 i • 

We expand j(y) in a Taylor series about x and substitute into the integral to obtain 

fa xx 2 [” a + ( 2 !<oc I 

i(y t) j( x ) ' i' ( x n 2 _ 2 1 2(2 x^rr^p 2 lOR L 1 f * Jr + hi ^ her order terms in a - 
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To first order we see 


y 


1 


1 f X 2 

2 \ a + x + 


log 


1 + 


a 


+ ( 2xa + a 2 
x 



We have used (6) to determine the evaluation point for all the intervals, namely 

(y. - yi-i) + (y , 2 ~yj-i) 1/2 j 


y. 


yw 


2 2\y - 2 - y 1 1/2 


log 


( 6 ) 


( 7 ) 


The required vector j. at the selected points given by (7) can be interpolated from an input 
tabular function representing the equatorial pitch angle distribution or the unidirectional flux per- 
pendicular to the field line, j x . If x m represents the mirror point value for a particle whose equa- 
torial pitch angle cosine is y , then Liouville's theorem states that 


i(y) - ii ( X J • 


From this it is seen that these two forms of input are handled in exactly the same way. 

The inverse matrix method of Reference 2 was not used for the L region below 3.0 to convert 
omnidirectional flux into unidirectional because the results with tabular inputs produced some 
bumps particularly where the spacing interval y i -y i _ 1 changed. In general this results because 
the slope of the input function is very important in the calculation. However, the continuity of this 
slope is not usually watched very closely when one is reading points off a graph, especially with a 
logarithmic scale. To understand the importance of the derivative, we have formulated the con- 
version in a different manner. The integral expression that relates J to j has been given by several 
authors (References 2, 3, and 4). We started with the expression given by Farley and Sanders 
(Reference 4). 


and using the equations 


_ R 3/2 r B 


) - 


f c dB d / J \ 

2.* J B (B-BJV 2 dB[Bj 


B o 

1 - x* 

Bp 

1 - y 2 

B o = B p 
1 - x c 2 1 - y c 2 
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have obtained 


i(y) 


1 r y c (i- x »)i/2 

2" 2 (i-y 2 )J y (x 2 -y 2 ) 1/2dx 


[(l-x 2 ) j] dx . 


The simpler expression given by Roberts (Reference 2) 


(8) 


j( y ) 2n 2 J x=y (x 2 -y 2 ) 1/2 dx x2 ) 1/2j ] dx 

can be integrated by parts to show that (8) and (9) are equivalent. 

By breaking up the x and y range into intervals, as has been discussed previously, one can 
write (8) as 





(l - x 2 ) 1/2 d 

( X 2_y i 2 i ) 1/2 


[(l-X 2 )j]dx . 


The assumption that j(x) can be approximated over the ith interval by 


Ji (x) 



c . e 

X 


-B 0 / b i (l -- 2) 


( 10 ) 


( 11 ) 


is known to be reasonable because all the data are plotted before input to the computer is made. 
Enough points are fed in to make interpolation in this form valid. Here c. and b. are constant 
over the interval. Differentiating (11) with respect to x and substituting into (10), we find 

2772 (i-Vk-i) i( y k-i) = 


N r * 

1. 


Jj (X) 


2x (l - x 2 ) 1/2 

(x 2 - y k -i ) 1/2 


2xB„ 


b s (l-x 2 )*/ 2 (x 2 - y ) 


V 2 


dx 


( 12 ) 


Again if the intervals are small enough, J 4 (x) can be removed from the integral. We have defined 
the D and E matrices by 


'ik 


2772 ( 


1 r x i 

1 -yu-i) Jx. 


2x(l - x 2 ) 1/2 

._ 1 (x 2 -y k 2 i ) 1/2 


dx 


i > k 




i < k 
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Bo p 

(! “Vk-i) J Vl 


2x dx 


/ 27T 2 (1-Vk-i) Jx.., (l ~ x 2 ) 1/2 ( x 2 ~ y ^. 1 ) 1/2 ' 


i > k 


i < k 


and the vectors 


F i B Ji( x i) 


F;' S 


Ji 

b. “ dB 


Jk " j( y k-i) • 


If these are used, then (12) can be written as the matrix equation 


N 

i k = £ (DikFi+EikFi ') 

i=0 


Using 


X 2 _ 2 

i y k-l 

l'X 2 


X 2 2 

i-1 y k-l 

1 - X 2 
1 x i-l 


1/2 


1/2 


the matrices can be expressed as 


1 ( p 


2n 2 \p 2 + 1 


+ tan 1 p - 


q 2 + 1 


tan 1 q ; P 2 > 0 


P 2 < 0 


< 


772 1 1 -yk-i) 


(tan 1 p - tan” 1 q ) ; p 2 > 0 


-0 ; p 2 < 0 ■ 

The two matrices were constructed as 45 by 45 triangular ones at each of the L values of interest. 
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The conversion of flux from one form into the other requires a knowledge of the omnidirec- 
tional flux at all points along the field line or of a complete equatorial pitch angle distribution. 
Most measurements only supply a portion of the required range. The most common case is shown 
in Figure 1. In order to make the conversion of these data, extrapolated data shown by dotted lines 
are also fed into the computer. The computer is made to calculate separately for the real data 
and extrapolated data. If the converted flux is composed of more than 75% real data, it is treated 
as real data in the new form; otherwise it is treated as extrapolated data. 

All of the flux data below L = 3.0 have been converted by this program so that J(B, L) and 
j 1 (B, L) exist on punched cards in our standard card format which is explained in the next section. 

For arbitrary L values we have employed the method given by Roberts (Reference 2) except 
that we have made a variable- spaced 45 by 45 matrix 


a 


ik 



i > k 


i < k 


where 


a i y c - yi 
a k y c = x k 

instead of a uniformly spaced one. The omnidirectional flux is obtained from 

i=l 

and the unidirectional flux is obtained from 

j ( y x-0 = t : ik lG (*i) 

i- 1 


where a i 1 is the inverse matrix of a , . 

ik ik 


PUNCHED CARD FORMAT 


We will describe the format we have chosen to store the trapped radiation data on punched 
cards since others may wish to use these cards. A given set of cards containing the data function 
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F(B, L) is comprised of a partitioned subset of cards; each subset corresponds to a fixed L value. 

At the head of each subset there is a header or title card followed by the data cards for that L value. 
Thus an entire deck of cards covering all L values is ordered so that the first card is a header 
card which is followed by its data cards; then the same sequence is repeated at each L value. The 
formats of the header card and data card are given below. 


Header Card Format 


Columns 1-60 Used to describe the data in alpha numeric symbols. 


Columns 61 - 66 Not used at present. 

Column 67 A zero (0) indicates that the data does not cover the whole range of B values. 

A one (1) indicates a complete set of data which spans from the equator to 
the atmospheric cutoff. 

Columns 68 - 69 Give the number of data cards behind the header card. 


Column 70 Not used. 

Column 71 Denotes the type of data with the following code: 

0 - Extrapolated data 

1 - Omnidirectional flux 

2 - Unidirectional flux 

3 - Unidirectional flux perpendicular to field line, j L 

4 - Counting rate, omnidirectional 

5 - Counting rate, unidirectional 

6 - Counting rate, perpendicular 

7 - Composite environment, omnidirectional flux 

8 - Composite environment, unidirectional flux 

9 - Other 

The composite environment refers to those produced in this program such 
as AE1, API, etc. 

Column 72 Denotes the particle and approximate energy threshold with the following code: 

0 - Protons <0.1 Me V 

1 - Protons 0.1 - 1.0 MeV 

2 - Protons 1-4 MeV 

3 - Protons 4-30 MeV 

4 - Protons >30 MeV 

5 - Electrons <30 keV 

6 - Electrons 30 - 300 keV 

7 - Electrons 0.3 - 1.0 MeV 

8 - Electrons >1.0 MeV 

9 - Unspecified at present 
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Columns 73 - 75 Identification of the source of the data. The first digit gives the experimental 
group or groups and the last two digits will index their data. The code for 
the experimental groups is: 

000 - 099 - Aerospace Corporation 

100 - 199 - Air Force Cambridge Research Laboratory 

200 - 299 - Bell Telephone Laboratories 

300 - 399 - Goddard Space Flight Center 

400 - 499 - Lockheed Missile and Space Company 

500 - 599 - University of California at San Diego 

600 - 699 - State University of Iowa 

700 - 749 - Applied Physics Laboratory, Johns Hopkins University 

750 - 799 - Rice University 

800 - 899 - Composite environment 

900 - 949 - University of Minnesota 

950 - 999 - University of California (including Lawrence Radiation Laboratory) 
Of course as other groups begin to contribute data, new assignments of num- 
bers will be made. 

Columns 76 - 79 Give the L value as x. x. 

Column 80 Not used. 

An example of a header card is shown in Figure 2. 

Data Card Format 

Columns 1-60 Contain four pairs of values of the form (B, F) for a fixed L value. Each pair 
has a field of fifteen columns; the first six columns are for B and the last 
nine are for F . The format B is xxxxx and the format for F is either xxxEyy 
or . xxxE+yy. The exponent is understood to be positive when it is not punched. 

Columns 61-70 Not used at present. 

Columns 71-78 Are same as header card. 

Columns 79 - 80 Used for indexing the data cards behind the header card. 

Two examples of data cards are shown in Figure 2. 

ENERGY SPECTRUM 

We have found it useful to represent the differential energy omnidirectional flux, J(E, B, L), 
as a product of two functions 

J(E, B, L) = F(B, L) 
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We call F(B, L) the distribution function and the term in the brackets is the differential spectrum 
normalized such that the integral spectrum 


N( >E, B, L) 


- f 

J F. 


N(E'. B, L) dE' 


r „ ( e’ . b, 

J E 

Je, 


L) dE' 


L)dE' 


is equal to 1 for E - Ej . We will suppress the B and L variables and write N(E, E x ) for the differ- 
ential spectrum andN(>E, E x ) for the integral spectrum normalized in this manner. Using (13) and 
this notation, the integral energy omnidirectional flux is written as 


J(>E, B, L) = F(B, L) N (>E, E x ) 


and we see that the distribution function represents the flux above energy E x . 

Because measurements are made with instruments of different spectral characteristics, it is 
imperative to have some knowledge of the energy spectrum in order to compare different sets of 
data. One is either provided with some direct spectral measurements or can obtain this informa- 
tion by an iterative process with several sets of data if the spectral characteristics of the various 
instruments are known. To illustrate the importance of this spectral information in making model 
environments, we will discuss the common threshold detector which has supplied much of the trapped 
radiation data. 

The efficiency versus energy curves for two of the electron detectors flown on Relay I are 
given in Figure 3 (References 5 and 6). Although these two detectors were actually directional 
ones, we will proceed with the discussion in omnidirectional form. The observational quantity of 
such a device is the counting rate, c , which is converted into a flux by the relationship 


F(B, L) 


c(B, L) 


/•CD 

e(E)N(E, E,) 
J 0 


dE 


Here g is the geometrical factor and we have expressed all of the energy dependence of the detector 
in the e(E) term. In order to understand more fully the limitations of the threshold detector, we 
consider the idealized efficiency curve given in Figure 4 and assume an exponential spectrum 


n(e,e t ) = 


exp 



E 


o 
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vjhich is normalized in our usual way at E T . It is easy to show that average efficiency 


e 


f 

J n 


e(E) N(E, E t ) dE 


% A 

€ o ~K Slnh e7 


in this case. We see that 


e 




1 . 


For this example 7 is constant to second order and the detector is insensitive to spectral changes. 
It is interesting to note that had we chosen to call the threshold energy of this idealized detector 
E t ± A, then the average efficiency would be 


±^0 


A sinh E„ 


or for A/E 0 « 1 , 


e = e 0 (l±VE 0 ) 


which is only constant to first order. Thus it is better to associate the threshold energy of a 
detector with the approximate mid point of the rising efficiency curve when the spectrum is not 
accurately known. 

If one uses a power law spectrum 


N(E, E x ) 


' E \~< b+2 > 


(b+ 1) . 
e t 


b > 0 


in the preceding analysis, then 


e 



/b + n 


1 1 1 

\ b ) 

_(1 

-Y) b (i + y) b J 


1 

1 1 , i 


27 

_(l-y) b (l+y) b+1 J ‘ ( l + y) b+1 
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where y - V E T . If y « l, then 


e 



1 + 


(b + 1) (b + 2) y 2 
6 


It is apparent that e(E) must be known for the different measurements in order to obtain a 
spectrum or to compare the data in detail; unfortunately this information is not always available. 


DISTRIBUTION FUNCTION 

When a given model spectrum is obtained (by whatever means are available), the various 
measurements can be converted to some common energy value. The energy value chosen is 
usually that of some extensive measurement or near several different measurements so that 
errors in the model spectrum will not produce large errors in converting fluxes to the common 
energy point. If an experimental measurement of flux above energy E 2 is denoted by 

J A >E r B ' L ) 


the quantity 


F m (B, L) 


J m (>e 2 . b,l) 

N(>E 2 , E t ) 


represents the data converted to a flux above energy E x by the model spectrum. If a measurement 
covers an energy band E 2 -E 3 , then the conversion is made by 


(B, L) 


J m (e 2 -e 3 , b.l) 
n(>e 2 , e,) -n(>e 3 , e,) • 


Following this procedure with all the available measurements, one obtains a series of curves at 
each L value. Some type of mean value is taken from this and becomes the F(B , L) for the model 
environment. It has been found necessary to draw the B-L flux map presentation in order to draw 
curves through the scatter of points at each L value in a somewhat consistent fashion. Examples 
of distribution functions and the converted data used in making them will be shown later in connec- 
tion with each model environment. 


ORBITAL INTEGRATIONS 


Having the environment in terms of the spectral and distribution functions is very useful for 
comparing with new experimental data or for designing space experiments. One of the most useful 
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wpys of displaying the data for engineering and planning purposes is by giving the average ac- 
cumulated flux that a satellite will receive in orbiting the earth. This is generally called the 
orbital integrated flux and will be denoted by 


«>e> = f T L (>E 1 L)d ’ • 

J o 

The calculation is performed on a computer using a series of standard subroutines. Many of these 
have been described elsewhere (Reference 7) so a detailed discussion will not be given here. An 
orbit program calculates the altitude, latitude, and longitude (h, 0,4>) at various times along the 
orbit. B(h, 0 , <p) is calculated using the 1960 Jensen and Cain coefficients (Reference 8), and then 
L is calculated. The flux for thisB, L point is then obtained and multiplied by the time interval 
between successive points. Proceeding in this fashion and accumulating the flux-time interval 
product, 4> is obtained. The time intervals used must be small enough so that the integral is 
adequately done by the trapezoidal rule and the total time, T , should be long enough to obtain the 
desired average. 

Since the manner in which the flux map is interpolated for a given B, L point varies, we will 
describe how it has been done for the orbital integrations presented later in this work. Flux values 
are stored for a series of B values at discrete L values. We denote the L values by L i and order 
them so that L. <L i+1 ; the B values, B u , are ordered so that B u < B iJ+1 . The B values are not 
necessarily the same on each L shell, i.e., f B i+1 s . To obtain the flux at any point B', L' , a 
four- point interpolation on the logarithm of the flux is performed. 

The flux at the point (B' , L' ) , where 

L, 1 L' < L i+I> 

*ij < B' < B i>j+lj 


and 


B i+1>k £ B' < B i+lik+1 


is obtained by 


F(B', L') = exp^(B'.L') . 


Here 




*(B',L H 1 )-l(B', Li ) 


L i+1 ~ L i 


(L'-Li) 
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*(B'. Li) = ^(Bii-Lj + 


^( B i,j+1’ L i)"^(B ij .L i ) 






and 


^(B » L i+l) ^(Bi+l.k* L i+l) 


+ 


^( B i+l,k+l’ L i+l) ^( B i+l,k’ L i+l) 


®i + l,k+l ®i+l f k 


(B' -B, 


We have used a standard set of circular orbits between 150 and 6000 n mi altitude to do the 
orbital integrations for the inner zone environments. For altitudes below 4000 n mi, a time in- 
terval of one minute was found to be adequate and a total time of one day will average out longitudi- 
nal asymmetries produced by the anomalies in the magnetic field. For the altitudes 4000 through 
6000 n mi, a two- minute step size and T = 2 days were used. These 76 orbits which comprise the 
standard set have been retained on punched cards so that calculations of b and L will not have to be 
repeated when integrations are done with new environments. For the low altitude orbits there are 
geographic regions where the flux is known to be negligible; a geographic test is made to determine 
if the flux should be set equal to zero without calculating B , L. This procedure reduces computer 
time considerably for low altitude orbits and has been used independently by others (Reference 9). 
An exclusion test for large L values is also used to reduce computer time when high inclination 
orbits are calculated. This is done by calculating L for a simple dipole and requiring thatL be 
less than some desired value to continue the calculation at that point. 


ELECTRON ENVIRONMENT AE1 

The first model environment produced in this program was the inner zone electron model 
designated AE1. This was based on data available to us by June 1964. The main measurements 
that were used in constructing AE1 are shown in Table 1. We shall use the code shown in the last 
column in referring to the data hereafter. Since it was desired to present an environment as up- 
to-date as possible, data taken during late 1962 that presented evidence of injected electrons from 
the Russian detonations and their subsequent decay into the natural background have not been used. 

At the time this environment was made, only the efficiency curve for E4 was available. The 
decay of the Starfish electrons during the time period of the various data also made it difficult to 
obtain good spectral information from the various threshold detectors. The model spectrum was 
obtained mainly from the detailed spectral measurements of E5 and E7. E7 covered the 1.17 to 
1.70 L region, whereas E5 covered both the inner and outer zones in late 1962. AtL ^1.25 and 
1.7, the two measurements E5 and E7 were in good agreement. Above 2.0 MeV at all L values, an 
exponential spectrum gave a reasonable fit. The fact that E5 did not show a very large B dependence 
for the energy spectrum and the different energy threshold detectors exhibited similar shapes 
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Data Used in Making AE1 Environment. 
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along a given L shell resulted in our using a model spectrum that depended only on L. An inde- 
pendent examination of E5 yielded similar results (Reference 9). In the region 1.25 <L< 1.7 
and E < 2 MeV, E5 and E7 were not in agreement. The reason for the differences are not clear but 
may be due to the energy dependence of Starfish decay (Reference 16), natural electron changes, or 
a large B dependence in the region where E7 was obtained. A compromise between the two meas- 
urements was used below L = 1.7, and E5 was followed closely at the higher L values although the 
effects of the Russian injections may have been present. 

In order to examine the behavior of the high energy portion of the spectrum, we plotted the 
exponential parameter E 0 for some of the electron data (Figure 5). Since various curves could be 
drawn through the E5 data, we have presented the E 0 obtained from the curves given in Reference 
11 as well as our own. The 0.5 MeV (El) and > 5.0 MeV (E2) measurements were used to obtain an 
E 0 although this was not expected to agree quantitatively with the other values because the spectrum 
is not exponential over the whole energy range. One sees the high energy portion of the spectrum 
gets soft as L increases from 1.2 to 1.8 and then becomes hard in the slot region L = 2 to 3. 

Although it is not completely clear that this behavior in the slot continued through the summer of 
1963, other Relay I measurements of Brown et al. (Reference 6) indicated no spectral change out 
to L = 2.3 in this time period. The solid line in Figure 5 shows the L variation of E 0 that was 
chosen; the increase of E 0 above the 0.5 MeV minimum was delayed until L = 1.9 to obtain better 
agreement between E3 and E4. The differential spectrum obtained from these considerations is 
shown in Figures 6 through 9 and is given in Table 2 along with the integral spectrum. 

Using this spectrum, all the data have been converted to >0.5 MeV so that a distribution func- 
tion could be obtained. We show the results of this conversion in Figures 10 through 23 in the 
following manner. The individual measurements are shown as points at specific B values and the 
solid line represents the environment chosen. We show these at many L values because this is 
the best way to display the agreement of the data and to show those regions where no data existed. 
The E3 data were used initially with no energy correction but they have been plotted here using 
the efficiency curve shown in Figure 3. The e for both detectors using the model spectrum is 
given in Table 3 for a 0.5 MeV threshold energy. For E4, 1.0 and 1.7 MeV are also shown. It is 
seen that the higher energy detector shows larger changes to the L dependent spectrum the lower 
the threshold energy. All the other threshold data were converted to >0.5 MeV on the basis of the 
energies given in Table 1. The E5 data were not converted to omnidirectional flux, but it is shown 
where the two data overlap as the El data reduced by a factor of two are given in Reference 11. 

For L values above 2.5 only El and E2 were available, and El was used as being representative of 
the fluxes in this region although there are time changes that occur. 

An examination of all these plots shows that the data spread over a considerable range and 
it is not really possible to see a systematic decay. Clearly the threshold for E2 depends on the 
spectrum as West, et al. (Reference 11) have pointed out. The agreement of E3 and E4 is reasonable 
up to L = 1.9; however to bring them in agreement at larger L values would require a considerably 
softer spectrum. The fact that the Brown, et al. detectors showed no spectral change in this time 
period makes this region a little uncertain. More recent data suggest considerable change in the 
high energy portion of the spectrum in this region. 
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Energy Spectrum- AE1. 
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Table 2 (Continued) 
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Table 3 


Average Efficiency for E3 and E4 Detectors. 


L 

Detector E3 

Detector E4 

6 0-5 

6 0-5 

j* 1 
6 

6 1.7 

1.2 

0.93 

6.28 x 10 ~ 2 

7.95 x 10' 2 

12.8 x 10' 2 

1.3 

1.00 

5.39 x 10 -2 

7.41 x 10 - 2 

13.1 x 10' 2 

1.4 

1.01 

4.35 x lO- 2 

6.68 x 10-2 

13.4 x 10 2 

1.5 

1.01 

3.51 x 10' 2 

5.94 x 10- 2 

14.1 x 10' 2 

1.6 

1.09 

2.72 x lO- 2 

5.89 x 10-2 

14.3 x 10 2 

1.7 

1.23 

1.50 x 10 '2 

4.93 x 10-2 

17.2 x 10' 2 

1.8 

1.44 

1.28 x 10'2 

5.52 x 10" 2 

15.4 x 10' 2 

1.9 

1.45 

1.27 x 10 "2 

5.92 x 10" 2 

14.5 x 10' 2 

2.0 

1.29 

2.56 x 10-2 

8.29 x 10" 2 

16.7 x 10' 2 

2.1 

1.21 

4.13 x 10"2 

9.18 x 10' 2 

14.6 x 10' 2 

2.2 

1.16 

5.50 x 10-2 

8.90 x 10' 2 

12.0 x 10' 2 

2.3 - 2.7 

1.03 

8.65 x IQ -2 

9.46 x 10" 2 

11.7 x 10 ' 2 

2.7 - 3.0 

1.07 





We tended to normalize the environment near E3 because this was a later time measurement 
and fell in the middle of much of the converted data. At higher L values the model falls closer to 
El. At low altitudes electrons are wiped off in the South Atlantic anomaly as they drift eastward 
and the depleted regions are repopulated during the remainder of the drift around the earth (Refer- 
ences 13 ; 17, and 18). This produces a longitudinal variation that becomes pronounced for mirror 
points whose minimum altitudes lie deep in the atmosphere. The low altitude fluxes also vary with 
magnetic activity based on more recent data (Reference 18). In this first electron environment we 
have not attempted to represent these effects in detail and have bounded the map at low altitudes 
by making the lowest flux contour (10 3 ) follow the lineh min = 0 km over much of the B-L space. 
Clearly there are both time and longitudinal variations in this region. 

The final F-B curves were obtained by comparing with the B-L flux map representation and mak- 
ing adjustments consistent with the data until smooth curves were produced in both representations. 
The final results are shown in Figures 24 through 26 as F-B, B-L and R-A. flux maps. The numerical 
values are available on punched cards and have been distributed to many organizations. 

The orbital integrations were performed with this environment by the method described in 
Section VII. The spectrum was handled in the following way. At each L value 


Lj = 1.1 + O.li 
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up to L = 3.0, the fraction of the particles in each of eight energy bands as given in Table 4 was 
specified. The computer used this as the spectrum over the whole interval 


h i ' 0.05£L<L i + 0.05 . 


The final results are given in Figure 27 forE > 0.5 MeV. The results for all the energy bands have 
been distributed in tabular form. The equatorial pitch angle distribution or j obtained by the 
conversion program described in Section III is shown in Figure 28. 


Table 4 

Fraction of Electrons in Energy Bands (MeV) - Map AE1. 


1.25 < L < 1.35 


1.35 < L < 1.45 


1.45 < L < 1.55 


E(I) 

E ( I +1) 

0.50 

1.00 

1.00 

2.00 

2.00 

3.00 

3.00 

4.00 

4.00 

5.00 

5.00 

6.00 

6.00 

7.00 

7.00 

00 

E ( I ) 

E ( I +1) 

0.50 

1.00 

1.00 

2.00 

2.00 

3.00 

3.00 

4.00 

4.00 

5.00 

5.00 

6.00 

6.00 

7.00 

7.00 

00 

E ( I ) 

E ( I +1) 

0.50 

1.00 

1.00 

2.00 

2.00 

3.00 

3.00 

4.00 

4.00 

5.00 

5.00 

6.00 

6.00 

7.00 

7.00 

00 


L < 1.25 

0.2100E-00 

0.4000E-00 

0.2120E-00 

0.9620E-01 

0.4410E-01 

0.2070E-01 

0.9200E-02 

0.7800E-02 

1.65 < L < 1.75 

0.6960E-00 

0.2550E-00 

0.4240E-01 

0.5710E-02 

0.7680E-03 

0.1050E-03 

0.1470E-04 

0.2300E-05 

2.15 < L < 2.25 

0.3820E-00 

0.2500E-00 

0.2000E-00 

0.9050E-01 

0.4210E-01 

0.1910E-01 

0.8850E-02 

0.7450E-02 


0.2740E-00 

0.4210E-00 

0.1920E-00 

0.7140E-01 

0.2610E-01 

0.9700E-02 

0.3700E-02 

0.2100E-02 

1.75 < L < 1.85 

0.7680E-00 

0.1850E-00 

0.4070E-01 

0.5450E-02 

0.7340E-03 

0.1000E-03 

0.1380E-04 

0.2200E-05 

2.25 < L < 2.75 

0.8800E-01 

0.3140E-00 

0.3220E-00 

0.1480E-00 

0.6830E-01 

0.3230E-01 

0.1460E-01 

0.1280E-01 


0.3490E-00 

0.4290E-00 

0.1620E-00 

0.4360E-01 

0.1190E-01 

0.3270E-02 

0.8950E-03 

0.3350E-03 

1.85 < L < 1.95 

0.7860E-00 

0.1640E-00 

0.4320E-01 

0.5870E-02 

0.8060E-03 

0.1070E-03 

0.1460E-04 

0.2400E-05 

2.75 < L < 3.05 

0.1100E-00 

0.3560E-00 

0.2800E-00 

0.1320E-00 

0.6270E-01 

0.3120E-01 

0.1450E-01 

0.1360E-01 


0.4090E-00 

0.4440E-00 

0.1250E-00 

0.1880E-01 

0.2740E-02 

0.3920E-03 

0.6140E-04 

0.6600E-05 

1.95 < L <2.05 

0.6910E-00 

0.2040E-00 

0.7720E-01 

0.2050E-01 

0.5390E-02 

0.1410E-02 

0.3680E-03 

0.1320E-03 


1.55 < L < 1.65 

0.5380E-00 
0.3510E-00 
0.9580E-01 
0. 1320E-01 
0.1720E-02 
0.2440E-03 
0.3110E-04 
0.4900E-05 

2.05 < L < 2.15 

0.5500E-00 

0.2370E-00 

0.1360E-00 

0.4840E-01 

0.1810E-01 

0.6600E-02 

0.2450E-02 

0.1450E-02 


UPDATING OF THE ELECTRON ENVIRONMENT 

An updating of the inner zone electron environment will be undertaken in the near future as 
soon as additional information is available. There will be spectral data from the EGO and P-11 
satellites as well as extensive mapping from Relay II, Injun IV, and Explorer 26 in the near future. 
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Some recent work by Walt (Reference 19) and Walt and McDonald (Reference 20) has shown that 
atmospheric scattering and energy loss explain the decay of the Starfish residue very well for 
L < 1.25, but theory and observation depart radically at higher Lvalues. Experimental decay life- 
times for ^2 MeV electrons have been reported by Van Allen (Reference 21) and for >1.2 MeV by 
Bostrom and Williams (Reference 16). We show these results in Figure 29. Hess (Reference 22) 
has used the lifetimes quoted in Reference 16 to decay his E8 grid from November 1962 to January 
1965, and has cautioned correctly that no account has been taken of the energy dependence of the 
decay nor is one certain when the other energies will reach the natural background. 

We have used the AE1 spectrum to convert 
the Bostrom and Williams >1.2 MeV measure- 
ments at L = 1.2, 1.3, 1.4, 1.5, 1.6, 1.7 to >0.5 
MeV using the efficiency curve given in Refer- 
ence 23. Examples of these results are given 
in Figures 13 and 14. With these results and 
the measured lifetimes we have calculated the 
time between their later measurements and 
AE1; the resulting time intervals are given in 
Table 5. It is seen that the ratio of their data 
and AE1 vary with L. The average interval over 
all L’s would place the time of AE1 around 1 
March 1963. Based on the normalization of AE1 near the E3 measurement, we would expect the 
date of AE1 to be 1 July 1963 and recommend using this date for decay purposes. The spread be- 
tween the API data and AE1 after correction for decay is about the same as much of the data in 
Table 1. Using the Bostrom and Williams decay values, we have performed orbital integrations 
for AE1 decayed in July 1, 1965. The results are given in Figure 32 along with the original results 
from AE1 for comparison. 


Table 5 


Decay Interval: AE1 Map and Bostrom- Williams 
Data Days 64-70, 1964. 


L 

Time Interval (Days) 

1.2 

262 

1.3 

347 

1.4 

402 

1.5 

418 

1.6 

414 

Avg. 

369 


THE INNER ZONE PROTON ENVIRONMENT 

The inner zone proton flux is relatively stable in time although some changes at low altitudes 
occur over the solar cycle because of atmospheric changes. The magnitude of this effect seems to 
be less than that predicted by Blanchard and Hess (Reference 24) probably because of errors in the 
atmospheric model chosen. McHwain (Reference 25) has observed a radical decrease in the high 
energy component >34 MeV for L > 2.0 following the large magnetic storm of 23 September 1963 
and is observing the slow return of this to prestorm values (Reference 26). Farther out in the 
magnetosphere, the proton distributions are more easily affected by magnetic disturbances but in 
general they are more stable than the electron fluxes. 


Although emulsions have been used to make some spectral measurements at low altitudes, 
most of the information concerning the distribution and energy spectrum has been obtained from 
threshold detectors which usually measure within an energy window. The model environments 
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described in this section will be for energies above 4 MeV. The lower energy protons which extend 
throughout the outer zone will be treated at a later date. The satellite data that have been used in 
constructing the proton maps are given in Tables 6 and 7. We shall refer to the data by the code 
name given in these tables. 

For certain energy intervals and for large regions of space the integral spectrum can be 
represented by an exponential spectrum. To illustrate this point we show the Naugle and Kniffen 
data (Reference 27) at high B values in Figure 33 and the equatorial spectrum obtained on Relay I 
(Reference 28) in Figure 34. These both suggest that the sum of two exponentials would be good 
over much of the inner zone. We have attempted to represent the environment in this fashion but 
have found this is too restrictive to describe the whole spatial and energy domain we desire to 
include. Instead, we have broken up the environment into four energy bands: 4 to 15, 15 to 30, 

30 to 50, and >50 MeV and have used an exponential spectrum within each band. The energy param- 
eter, E 0 , is a function of both B and L so we are representing the omnidirectional integral flux as 


J(>E, B, L) 


F(B, L) e 


E 1 -e/e 0 (B,L) 


By selecting a known energy function, we can still represent the environment in terms of functions 
of two variables rather than the more cumbersome general function of three variables. 


Proton Map API 

We will discuss the various environments in the same sequence in which they were derived. 

The data above 30 MeV were studied first. We have used the published energy intervals or thresh- 
olds for each srt of data. There were not many efficiency curves available and even if there were, 
the effort was not considered necessary in this first attempt at describing the environment. We 
did use the energy interval 50 to 75 MeV for the Telstar data although it has been referred to as 
>50 MeV in some cases. This choice was based on the published detector efficiency in a typical 
direction (Reference 33) and this choice gave better agreement with other data in all regions ex- 
cept where the spectrum was very hard. 

The function E 0 (B, L) obtained for this region is given in Figure 35 and we show the data con- 
verted to flux >34 MeV in Figures 36 through 52. This conversion is described in Section IV. We 
have selected E t = 34 MeV because the P6 data are the most extensive at the higher L values 
where there are not enough data to obtain a good spectrum. 

At low L and high B values the 1958 to 1961 data (PI, P2, and P3) seem to be separated from 
the 1962 data (P4) and probably reflect the solar cycle changes that are predicted. An examination 
of Figures 37 through 42 shows the effect diminishing as L increases or, more appropriately, as 
h min increases. We have made the model conform with the later data because of the time period 
assigned to the environment. The divergence of the Telstar data at high B values where the spectrum 
is hard is probably due to an energy dependence of energy window values. Of course this is probably 
true of the other detectors also. At any rate, we have selected the lines shown in the figures as 
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being representative. The resulting distribution function is shown in the three representations 
commonly used to display environments in Figures 53 through 55. 

Proton Map AP2 

From the examination of Figures 33 and 34, one suspects that the spectral slope will change 
radically somewhere in the region between 6 and 20 MeV depending on the B, L coordinates. Since 
there were very little data in this energy interval, we decided to make the change at 15 MeV in- 
dependent of the spatial point. 

The AP2 environment was constructed by using P7 data and the API distribution function to 
obtain the E 0 function to be used in the 15 to 30 MeV interval. This function is shown in Figure 56, 
and comparison with Figure 37 shows that the two are very similar. In fact, it is possible to use 
AP2 over a wider energy range than 15 to 30 MeV without much error. The various forms of the 
distribution function are given in Figures 57 through 59. 


Proton Map AP3 

It is interesting to note that the spectrum gets softer on a given L shell as one approaches 
the equator, at least for L < 2.0 over the energy range 15 to 50 MeV. At the higher energies the 
opposite effect occurs on the basis of the data of Imhof, and Smith (Reference 34) and the recent 
measurements of Freden, Blake, and Paulikas (Reference 35). Therefore, an additional map had to 
be constructed to describe the high energy protons. 

There were very little data available above 60 MeV with which to obtain an environment. How- 
ever using the data of References 34 and 35 along with that of Gar mire (Reference 36) we have 
obtained E 0 for the high energy map which we have designated AP3. The distribution function for 
AP3 is produced by API for E> 50 MeV and the various representations are shown in Figures 60 
through 62. The E 0 is plotted in Figure 63. 

Proton Map AP4 

Although there are data available down to 1 MeV in the inner zone, the spectrum changes 
enough that a map AP4 only covering the 4 to 15 MeV energy range was made. Extensions to 
lower energy will be done in the near future. The data used in making AP4 are given in Table 7. 

An estimate of the proper spectrum around 4 MeV in the inner zone was made using the P12 - 
P14 data as well as the 1.1 to 14 MeV protons measured by Mcllwain, et al. (Reference 5) and the 
2.5 to 3.8 MeV protons measured by Brown, et al. (Reference 6). The differential spectrum near 
the equator is shown in Figure 63. It is clear than an exponential integral spectrum in this region 
is not completely correct, but it is adequate for our purposes as long as one considers E > 4 MeV. 
The comparison of this model with the data is shown in Figures 65 through 79. The discrepancies 
at the low L values are believed to be partly due to our choice of the spectrum. It is obvious from 
the figures that we have followed P13 closely; this was the only measurement available above 
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I* = 3.0 in the energy range. The three representations of the distribution function are shown in 
Figures 80 through 82 and the spectral parameter is plotted in Figure 83. The equatorial omni- 
directional fluxes calculated with the proton environment are given in Figure 84 for various 
energies. The perpendicular fluxes, j x , for the four maps obtained by the conversion program 
discussed in Section in are shown in Figures 85 through 88. 

Using these omnidirectional flux maps, the orbital integrations for the standard circular orbits 
were performed. The logarithmic interpolation described in Section VH was used for both F(B, L) 
andE 0 (B, L). Although the four maps API to 4 are most accurate over the energy bands described 
earlier, it is possible to obtain reasonable answers outside these bands. Since it is sometimes 
cumbersome to use more than one map for a specific purpose, we have computed the orbital in- 
tegrations for 23 energy intervals with each map. These intervals are given in Table 8; it is hoped 
that these results will be useful in selecting a map best suited for a given purpose. The results 
of these calculations are too extensive to be included here but will be available under separate 
cover. We provide an example of the printout in Table 9. We show some of the orbital integration 
results for each map in Figures 89 through 92. 


Table 8 


Energy Bands Used in Orbital Integrations. 


API 

AP2 

AP3 

AP4 

5-10 

5-10 

5-10 

1-2 

10 - 15 

10 - 15 

10 - 15 

2-3 

15 - 20 

15 - 20 ♦ , , 

15 - 20 

3-4 


_ _ Recommended 



20 - 25 

20 - 25 

20 - 25 

4-5 t 

25 - 30 

25 - 30 f** 0 

25 - 30 

5-6 1 

30 - 35 f 

30 - 35 

30 - 35 

6-7 

35 - 40 Recommended 

35 - 40 

35 - 40 

7-8 

40 - 45 Range 

40 - 45 

40 - 45 

8-9 

45 - 50 J 

45 - 50 

45 - 50 

9-10 Recommended 

50 - 60 

50 - 60 

50 - 60 

10 - 11 Range 

60 - 70 

60 - 70 

60 - 70 | 

11 - 12 

70 - 80 

70 - 80 

70 - 80 1 

12 - 13 

80 - 90 

80 - 90 

80 - 90 

13-14 j 

90 - 100 

90 - 100 

90 - 100 

14-15 | 

100 - 110 

100 - 110 

100 - 110 

15 - 20 

110 - 120 

110 - 120 

110 - 120 Recommended 

20 - 25 

120 - 130 

120 - 130 

120 - 130 Range 

25 - 30 

130 - 140 

130 - 140 

130 - 140 

30 - 35 

140 - 150 

140 - 150 

140 - 150 

35 - 40 

150 - 200 

150 - 200 

150 - 200 

40 - 45 

200 - 250 

200 - 250 

200 - 250 

45 - 50 

250 - 300 

250 - 300 

250 - 300 | 

50 - 60 

300 - oo 

300 - ® 

300 - co \ 

60 - 70 
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Table 9 

Example of Orbital Integrations Printout 
Orbital Integration Map API. 

Orbit Altitude: 150 n mi Total Time: 24 Hours Time Interval: 1 Minute 


Orbit 

Inclination 

Degrees 

0 

30 

60 

90 

E k 

E k + i 

*(0 

*(E ktI -E k ) 

<*( >E k ) 




*K) 

< K E k + r' E k) 

5 

10 

0. 

0. 

0.400E 06 

0.845E 04 

0.272E 07 

0.170E 07 

0.287E 07 

0.199E 07 

10 

15 

0. 

0. 

0.392E 06 

0.825E 04 

0.101E 07 

0.479E 06 

0.879E 06 

0.471E 06 

15 

20 

0. 

0. 

0.384E 06 

0.807E 04 

0.534E 06 

0.147E 06 

0.407E 06 

0.120E 06 

20 

25 

0. 

0. 

0.376E 06 

0.788E 04 

0.386E 06 

0.5 19E 05 

0.287E 06 

0.360E 05 

25 

30 

0. 

0. 

0.368E 06 

0.771E 04 

0.334E 06 

0.225E 05 

0.25 IE 06 

0.145E 05 

30 

35 

0. 

0. 

0.360E 06 

0.753E 04 

0.312E 06 

0.127E 05 

0.237E 06 

0.848E 04 

35 

40 

0. 

0. 

0.352E 06 

0.736E 04 

0.299E 06 

0.927E 04 

0.228E 06 

0.650E 04 

40 

45 

0. 

0. 

0.345E 06 

0.720E 04 

0.290E 06 

0.784E 04 

0.222E 06 

0.571E 04 

45 

50 

0. 

0. 

0.338E 06 

0.703E 04 

0.282E 06 

0.715E 04 

0.216E 06 

0.530E 04 

50 

60 

0. 

0. 

0.331E 06 

0.136E 05 

0.275E 06 

0.132E 05 

0.2 HE 06 

0.991E 04 

60 

70 

0. 

0. 

0.317E 06 

0.130E 05 

0.262E 06 

0.123E 05 

0.201E 06 

0.928E 04 

70 

80 

0. 

0. 

0.304E 06 

0.124E 05 

0.249E 06 

0.115E 05 

0.191E 06 

0.876E 04 

80 

90 

0. 

0. 

0.292E 06 

0.118E 05 

0.238E 06 

0.109E 05 

0.183E 06 

0.830E 04 

90 | 

100 

0. 

0. 

0.280E 06 

0.113E 05 

0.227E 06 

0.104E 05 

0.174E 06 

0.788E 04 

100 

110 

0. 

0. 

0.269E 06 

0.108E 05 

0.2 16E 06 

0.989E 04 

0.167E 06 

0.748E 04 

110 

120 

0. i 

0. 

0.258E 06 

0.103E 05 

0.206E 06 

0.940E 04 

0.159E 06 

0.710E 04 

120 

130 

0. 

0. 

0.247E 06 

0.993E 04 

0.197E 06 

0.893E 04 

0.152E 06 

0.675E 04 

130 ' 

140 

0. 

0. 

0.237E 06 

0.949E 04 

0.188E 06 

0.849E 04 

0.145E 06 

0.642E 04 

140 

150 

0. 

0. 

0.228E 06 

0.908E 04 

0.180E 06 

0.808E 04 

0.139E 06 

0.610E 04 

150 

200 | 

0. 

0. 

0.219E 06 

0.398E 05 

0.172E 06 

0.349E 05 

0.133E 06 

0.263E 05 

200 1 

250 ' 

0. 

0. 

0.179E 06 

0.320E 05 

0.137E 06 

0.273E 05 

0.106E 06 

0.206E 05 

250 

300 

0. 

0. 

0.147E 06 

0.258E 05 

0.109E 06 

0.215E 05 

0.860E 05 

0.163E 05 

300 ; 

00 

0. 

0. 

0.12 IE 06 

0.12 IE 06 

0.882E 05 

0.882E 05 

0.697E 05 

0.697E 05 
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OMNIDIRECTIONAL FLUX 
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Figure 1 —Illustration of real and extrapolated data 
used in the conversion program. By extending the 
data, on both sides if necessary, the conversion 
of the flux is more correct near the end points of 
the real data. 
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Figure 2— Examples of header and data cards used to store 
the trapped radiation data. The two datacards show dif- 
ferent positions of the flux values within the field. See 
text for explanation. 
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Figure 3— Counter efficiency versus energy for two 
detectors flown on the Relay I satellite. 


Figure 4 — idealized threshold detector efficiency curve 
versus energy. The quantity e 0 is the maximum effi- 
ciency reached at any energy. 



1.0 1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 


L, EARTH RADII 


Figure 5 — The exponential energy parameter versus 
L for electron data. The solid curve is the varia- 
tion of E 0 chosen for the AE1 model environment. 
This model spectrum is exponential above 2.0 MeV. 
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Figure 10— Comparison of electron map AE1 with 
experimental data at L = 1 .2. The measured data 
are converted to the omnidirectional flux above 
0.5 MeV by means of the model spectrum. The 
coded points are identified in Table 1. 
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Figure 12— Comparison of electron map AE1 with 
experimental data at L = 1 .4. 
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Figure 1 1 —Comparison of electron map AE1 with 
experimental data at L = 1 .3. 
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Figure 1 3— Comparison of electron map AE1 with exper 
mental data at L = 1 .5. 
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| Ure 14— Comparison of electron map AE1 with experi- 
mental data at L = 1 .6. 


Figure 1 5— Comparison of electron map AE1 with experi 
mental data at L = 1 .7. 
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Figure 16— Comparison of electron map AE1 with experi- 
mental data at L = 1 .8. 

Figure 1 7— Comparison of electron map AE1 with exper 
mental data at L = 1 .9. 
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Figure 18— Comparison of electron map AE1 with experi- 
mental data at L = 2.0. 
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Figure 1 9— Comparison of electron mapAEl with experi- 
mental data at L = 2.1 . 
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Figure 20— Comparison of electron mapAEl with experi- 
mental data at L = 2.2. 


Figure 21 —Comparison of electron mapAEl with experi- 
mental data at L = 2.3. 
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Figure 25— The B-L flux map for AE1, 
E > 0.5 MeV. 
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Figure 26— The R-\flux map for AEI, 
E > 0.5 MeV. 
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igure 27— Orbital integrations with AE1, E> 0.5 MeV 
for circular orbits. 
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Figure 29— Apparent mean lifetimes of energetic electrons 
versus L. The results of several satellites are combined 
to show the general agreement of the measurements and 
to demonstrate the shorter lifetime for L > 2.0. 
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Figure 28— The flux versus B plot for j ^ or j(B 0 /sin 2 a 0 ) 
of AE 1 . 
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Figure 30— Comparison of AE1 and the APL meas- 
urements from satellite 1963 - 38c at L = 1 .2 and 
1.3. The data of Bostrom and Williams (Ref. 16) 
for Oct. 15-22, 1963 and March 4-10, 1964 are 
converted to > 0.5 MeV using the model spectrum 
of AE1 . 
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Figure 31 —Comparison of AE1 and the API measure- 
ments from satellite 1963 - 38c at L = 1.5. The shapes 
of the API data and AE1 are lesssimilar at this L value 
than those shown in Fig. 30. 



Figure 32— Comparison of orbital integrations with AE1 
for July 1, 1963 and July 1, 1965. The 1965 results 
were obtained by using the Bostrom-Williams decay 
constants and a time period of two years from the 
original AE1 environment. The angles shown are the 
orbit inclination angles. 
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Figure 33— Integral proton spectra obtained with the 
NERV experiment. The solid lines through the points 
at L = 1.79, 1.72, and 1.64 are the sum of two ex- 
ponential spectra. 
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Figure 34 — Integral proton spectra obtained at the 
equator with Relay I. The solid lines drawn through 
the points are the sum of two exponential spectra. 


EXPONENTIAL ENERGY 
PARAMETER E 0 
FOR PROTON MAP API 


2.4 2.3 

.5 1 / 2.2 


R.7M 


2.82.9 2., '2.0 



1.9 1.7 


0.04 0.08 0.12 0.16 0.20 0.24 

B, gauss 

Figure 35— The energy spectral parameter, E 0 , for pro- 
ton map API. This function is to be used over the energy 
range 30 - 50 MeV. 
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Figure 46 — Comparison of proton map API with 
experimental data at L = 2.2. 
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Figure 47 — Comparison of proton map API with 
experimental data at L = 2.3 
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Figure 48— Comparison of proton map API with 
experimental data at L - 2.4 
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Figure 49 — Comparison of proton map API with 
experimental data at L = 2.5. 
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Figure 50— Comparison of proton map API with 
experimental data at L = 2.6. 
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Figure 51 — Comparison of proton map API with 
experimental data at L = 2.7. 
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Figure 52— Comparison of proton map API with 
experimental data at L = 2.8. 
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Figure 53— The flux versus B plot of API, E > 34 MeV. 
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Figure 55— The R-\ flux map for API. The contours are 
the omnidirectional flux above 34 MeV. 
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Figure 56— The energy spectral parameter, E o/ for proton 
map AP2. This function is to be used over the energy 
range 15-30 MeV. 
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Figure 54— The B-L flux map for API . The contours are 
the omnidirectional flux above 34 MeV. 



Figure 57— The flux versus B plot of AP2, E >15 MeV. 















EARTH RADII map AP3. This function is to be used for energies greater 

than 50 MeV, 


Figure 62— The R- A flux map for AP3. The contours are 
the omnidirectional flux above 50 MeV. 
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Figure 64— Differential proton spectrum near the equator 
for inner zone protons. This is estimated from Relay I and 
Explorer 15 measurements discussed in the text. 
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igure 65— Comparison of proton map AP4 with experi- 


mental data at L = 1 .3. 
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Figure 67— Comparison of proton map AP4 with ex- 
perimental data at L = 1 .5. 


Figure 68— Comparison of proton map AP4 with ex- 
perimental data at L = 1.6. 





Figure 69— Comparison of proton map AP4 with 
experimentai data at L = 1.7. 



Figure 71 —Comparison of proton map AP4 with 
experimental data at L = 1.9. 



Figure 70— Comparison of proton map AP4 with 
experimental data at L = 1.8. 



Figure 72— Comparison of proton map AP4 with 
experimental data at L = 2„0 



Figure 73— Comparison of proton map AP4 with 
experimental data at L = 2.1 . 


Figure 74— Comparison of proton map AP4 with 
experimental data at L = 2.2. 





Figure 75— Comparison of proton map AP4 with experi- 
mental data at L = 2.3 


Figure 76— Comparison of proton map AP4 with experi- 
mental data at L = 2.4. 




Figure 77— Comparison of proton map AP4 with experi- 
mental data at L = 2.5. 


Figure 78— Comparison of proton map AP4 with expen 
mental data at L = 2.6. 



Figure 79— Comparison of proton map AP4 with experi- 
mental data at L = 2.7. 
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Figure 89— Orbital integrations with API, 
E > 34 MeV. 
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Figure 90— Orbital integrations with AP2, 

E > 15 MeV. 
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Figure 91— Orbital integrations with AP3, 
E > 50 MeV. 
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Figure 92— Orbital integrations with AP4, 
E > 4 MeV. 
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